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What determines the mass of the most massive star 
cluster in a galaxy: statistics, physics or disruption? 

M. Gieles^ 
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Abstract In many different galactic environments the 
cluster initial mass function (CIMF) is well described 
by a power-law with index —2. This implies a linear 
relation between the mass of the most massive cluster 
(Mmax) and the number of clusters. Assuming a con- 
stant cluster formation rate and no disruption of the 
most massive clusters it also means that M,„ax increases 
linearly with age when determining Afmax hi logarith- 
mic age bins. We observe this increase in five out of 
the seven galaxies in our sample, suggesting that Mmax 
is determined by the size of the sample. It also means 
that massive clusters are very stable against disrup- 
tion, in disagreement with the mass independent dis- 
ruption (MID) model presented at this conference. For 
the clusters in M51 and the Antennae galaxies the size- 
of-sample prediction breaks down around 10® A/q, sug- 
gesting that this is a physical upper limit to the masses 
of star clusters in these galaxies. In this method there 
is a degeneracy between MID and a CIMF truncation. 
We show how the cluster luminosity function can serve 
as a tool to distinguish between the two. 
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2 (e . g. lElmegreen and Efremovl 19971: Zhang and Fal 



fl993 iBik et all l2003l : IWhitmorel booa Ide Grijs et all 
120031 ). It is often assumed that the CIMF is "universal" . 
But is that really true? In this contribution we inves- 
tigate what the consequences are of a universal power- 
law CIMF with index —2. This function has several 
appealing properties that we use to make predictions, 
in particular for the most massive cluster in a galaxy. 

Lets consider the "universal" scenario, in which all 
cluster masses are drawn from the same CIMF: 
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For a < 2 the total mass diverges at high masses, so 
4>{M) has to be truncated at some upper value, which 
can be interpreted limit above which no clusters 
can exist. We will refer to this limiting mass as Mup 
and to the most massive cluster observed, i.e. the most 
massive cluster actually formed, as Mi„ax- By choosing 
Map >> Mmax the CIMF behaves as if it has no trun- 
cation and the value of Mmax is the result of sampling 
statistics. 

Afmax depends on the constant A in Eq. [T] and can 
be found by solving 



1 The most massive star cluster from 
statistical considerations 

In many studies of young star clusters the power-law 
is an ever returning function that is used to describe 
almost everything, where the only thing that varies is 
the index. The cluster initial mass function (CIMF) 
is one of the famous power-laws, with its index being 
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results in 
A/max for a = 2. The 
a depends on the num- 
ber of clusters formed (N) which is proportional to A. 
The larger N, the closer the statistically probable Afmax 



relation between M^p and Af, 



will be to A/up. The relation between A/„ 
be found from 



and N can 
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For a constant Af,nin we see from Eq. 2] that Af„ 
scales with as 



(6) 



Assuming a constant cluster formation rate (CFR), 
diV/dt=constant, the number of clusters in logarith- 
mic age bins, A'^bin; increases linearly with age, since 
7Vbi„(t) oc diV/dln(t) = tdN/dt and so Nun{t) oc t. 
The same holds for diV/dlog(t) apart from an addi- 
tional constant In(lO). The refore Mma^ in loga rithmic 
age bins scales with age as ( Hunter et al. 20031 ) 



(7) 



2 Variations in the cluster formation rate and 
mass independent disruption 

The derivation of Eg . [7| was based on a c onstant CFR 



and no disruption. Hunter et al. ( 2003| ) already no 



ticed that a power-law relation for the CFR with age 
(diV/dt oc f, with rj negative for a CFR that was lower 
in the past) would change the relation for Af^ax to 



(8) 



In a similar way we can add the effect of the 
mass indepe ndent dissolution ( MID) model that wa s 
proposed by iFall et all |2005h . IChandar et all |2006h . 
Whitmore et all |2007t l and Chandar (in these proceed- 
ings). These authors argue that 90% of the clusters are 
destroyed each age dex which, combined with a con- 
stant CFR, results in dN/dt oc for mass limited 
cluster samples. For other MID percentages we can 
write dN/dt oc tr^ with < A < 1 for MID percent- 
ages between 0% and 90%. Adding MID to Eq. H we 
get 
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From Eq.[H]it is visible that there is a degeneracy be- 
tween MID and a CFR that has been increasing from 



the past to now {jj < 0). Because of the mass indepen- 
dent nature of this disruption model it is impossible 
to distinguish between an increasing CFR and MID. 
For the u niversal 90% MID scen ario, which was pro- 
posed by Iwhitmore et all ( 2007 ). the predicted slope 
in a log(Afmax) vs. log(age) diagram is for all values 
of a in all galaxies. 

If A/max = A/up, i.e. a physical maximum to the 
cluster masses that is reached at all ages, Afmax is also 
constant with age. Hence there is also a degeneracy 
between the 90% MID model and a truncation of the 
CIMF. In §[l]we will show that the luminosity function 
is different in these two situations. In the following sec- 
tion we will compare empirical cluster ages and masses 
to the predictions for Afmax(i)- 



3 Comparison to observations 

3.1 Description of the data used 

From literature we collect ages and masses of star 
clusters in seven different galaxies: the Milky Way 
(solar neighbourhood), SMC, LMC, M33, M83, M51 
and the Antennae galax ie s. The data were taken 



Hunter et all (120031) . Mor a et al 



from iLamers et al.l (|2005[ ): iKharchenko et al 



( ioosh . 

(2007 , in p rep), 



Bastian et al.l (|2005r ) a nd IZhang and Falll (119991) re- 
spectively. We refer to ICieles and BastianI ( 2008h for 
more details on the data. 

3.2 Observed trends of log(Af,„ax) "vs. log(age) 

In Fig. [T] we show the trends of log(Arinax) with log(age) 
in bins of 0.5 dex as open circles. We also show the third 
most massive cluster in each age bin (A/max,3rd) as filled 
circles, which is expected to follow the same relation as 



ATrr 



but with less scatter. The full lines are fits of 



log(Afi„ax,3id) with log(age) over the first 100 Myr and 
the dashed lines consider an age range of 3 Gyr. 

For the solar neighbourhood, SMC, LMC, M33 and 
M83 the observed increase of log(Afi„ax,3rd) in the first 
100 Myr is consistent with the size-of-sample prediction 
without disruption (E g. [71). This rules o ut the long term 
90% MID proposed bv iwhitmore et al\ kooj ) as a uni- 
versal cl/uste7_di^^ for these galaxies. 
If the Whitmore et al\ \200'i ) model was correct, Mmax 
had to be constant with log (age) for all galaxies. Fig.\^ 
shows that this is clearly not the case. 

The trends for the interacting galaxies M51 and the 
Antennae are clearly different from the other five galax- 
ies. The flat slope beyond 100 Myr in M51 was inter- 
preted as a truncation of the CIMF around A/up ~ 
5 X IQ^Mq ICieles et al.ll2006al |bh. We note that M83 
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Fig. 1 Evolution of log(Afmax) in log(age) bins for the seven galaxies in our sample (open circles). The third most massive 
clusters (A/max, 3rd) in cach bin are shown as filled circles. Fits to log(Afmax,3rd) vs. log(age) on two different age ranges are 
shown as dashed and full lines. 



was only partially imaged so things could be different 
for this galaxy when all clusters are included since its 
CFR is probably similar to that of M51. 

The flat relation for the Antennae galaxies is consis- 
tent with 90% MID each age dex during a Gyr. Note, 
however, that if we attribute this result entirely to MID, 
i.e. we assume that the CFR has been roughly con- 
stant and that t here is no truncation o f the CIMF (the 
interpretation of lWhitmore et al.ll2007l ). then based on 
sizc-of-samplc effects the Antennae galaxies should have 
formed clusters with masses up to ~ 10^ Mq, but they 
have been destroyed by MID. This seems unlikely given 
what we know about other galactic mergers (i.e. that 
the major burst of star formation in the Antennae will 
happen as the nuclei coalesce, e.g. ICox et al. I l2006h . 



We try to put the results of Fig. [T] in a unified Mmax 
vs. N diagram, similar to wha t was done for the bright- 
est c l uster luminosity vs . N (Larsen 20021 : Whitmor3 



3.3 A universal cluster initial mass function? 

The galaxies in Fig. [T] are (roughly) ordered by increas- 
ing star/cluster formation rate. The slope of M^^x 
with age is getting smaller going from the solar neigh- 
bourhood to the Antennae galaxies. If cluster masses 
are drawn from a universal CIMF with a truncation 
at ~ 10^ Mq it is expected that for galaxies with a 
low CFR the CIMF is not sampled up to this mass 
(-^max << My^p) and thus Mmax follows from the size- 
of-sample prediction of Eq. [T] The truncation should 
then become more noticeable at high ages (larger age 
ranges) or in galaxies with a high CFR such as MSI 
and the Antennae galaxies. 



2003 ; Gieles et al. 2006a[ ) . Doing this for masses rather 



than luminosities is more difficult since we need to de- 
fine N as the number of clusters above a certain mass 
(Merit), with Afcrit tlic samc for all galaxies and all age 
bins. However, M^ax in the nearby galaxies is of the 
samc order as the minimum cluster mass we have avail- 
able for M5I and the Antennae and the minimum ob- 
servable cluster mass increa ses with age in sarn ples that 
are luminosity limited (e.g. iGieles et al.ll2007( ). 

From the available ages and masses we count the 
number of clusters with ages younger than 10® yr and 
masses above Merit, with log(Mcrit/MQ) = 2 for the 
solar neighbourhood, log(Mcrit/MQ) = 3 for the SMC 
and LMC, log(Mcrit/Mo) = 3.5 for M33 and M83, 
log(A'/crit/Mo) = 4 for M51 and log(A//crit/Af0) = 
5 for the Antennae. We then scale this number, 
N{M > Merit; i < 10® yr), to a formation rate of clus- 
ters with masses above 10^ Mq, diV(M > 10^ AfQ)/dt, 
by dividing A^(M > M^it; i < 10® yr) by 10® yr and 
multiplying it by Afcrit/10'^. The scaling with Merit 
is justified by Eq. El 

We can now multiply the width of each bin of each 
galaxy in Fig. [1] by the dN{M > 10^MQ)/dt of that 
galaxy to get an estimate of the number of clusters with 
masses above 10'^ Mq that has formed in each bin, again 
assuming a constant CFR. This allows us to plot all 
galaxies from Fig. [T] in the same figure. 
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Fig. 2 Mmax VS. Af(M > 10^ Mq) for all galaxies from 
Fig. [U The Af(M > 10^ Mq) values for each bin of Fig. [T] 
were obtained by multiplying the bin widths (At) from 
Fig. [T]by dN{M > lO^M0)/dt (see the text for details). 
The Mmax values are the same as in Fig. [1] for all galaxies. 
The dotted line represents a prediction for Mmax follow- 
ing from a power-law CIMF without a truncation (Eq. [5]). 
The dashed and full lines show a prediction based on a 
CIMF with Mup = 10*^ Mq and a Schechter type CIMF 
with M« = 10*^ Mq , respectively. 

The result is shown in Fig. [5] For Mmax smaller 
than ^lO^Mo (and N{M > 10^ Mq) < 10^) all 
galaxies show the same linear increase of Mmax with 
N{M > 10^ Mq), similar to that found for the brightest 
cluster luminosity. We plot the prediction for Afmax as a 
function of N from Eq.[5l i.e. for random sampling from 
a —2 power-law CIMF, as a dotted line. There is good 
agreement with the observations for Afmax ^ 10^ Mq. 
For higher Mmax the observed values flatten around 
IO^Mq. 

Using Eqs. [11^3] we solve for Mmax as a function 
of N{M > lO^Mg) for M^p = 10^ Mq (dashed line) 
and we do the same for a Schechter type CIMF with 
M* = 10^ Mq (full line). Both predictions including a 
truncation describe the data very well suggesting that 
the masses of star clusters follow the size-of-sample pre- 
diction only up to ~ 10^ Mq. Above that mass there 
are far fewer clusters observed than predicted by ran- 
dom sampling from a power-law CIMF without a trun- 
cation. 



4 The cluster luminosity function as an 
independent check 

Because the interpretation of the trends of log(Mmax) 
vs. log(age) could be affected by variations in the CFR 



and disruption (Eq. [9|) , we show how the integrated lu- 
minosity function (LF) will look like when the CIMF is 
truncated. When the mass function (MF) of a stellar 
population has a characteristic mass scale, then the LF 
of a multi-age population can not be inte rpreted as the 



MF due to the fading of cluster s with age (|Meurer et al 



199,4 iHogg and Phinnevll 19971 ) . However, if all clusters 
are drawn from a power-law distribution with index —2, 
without a truncation, the integrated LF should still be 
a power-law with index —2. Age dependent extinction, 
star-bursts or variations in the CFR and mass inde- 
pendent dissolution will only affect the normalisation 
constant, but not the power-law index. No matter how 
you normalise, a power-law plus a power-law is always 
a power-law. Mass dependent disruption will make the 
LF shallower, i.e. with an index larger than —2. 

However, the LFs of multi-age cluster populations 
seem to have power-law indices that are smaller (be- 
tween —2.5 and —2) (e .g. Dolphin and KennicuttI 2002 : 
Elmegreen et al.l l2002t iLarsenl l2002t ) with theTF be 



ing steeper at highe r luminosities (IZepf et al.l 11999 
Benedict etlH I2OO2I : iLarsenl l2002t IWhitmore et al 



1999 : lMengel et al.l2005HGieles et al.l2006aUbl : lHwang and Le^ 
2OO7I) . 

A slightly steeper LF results naturally from a trun- 
cated CIMF combined with evolutionary fading. If the 
CIMF is populated up to M^p at all ages, i.e. Mmax = 
A/up, then we have a constant number of the most mas- 
sive clusters at all ages: diV(Mup)/dt =constant. The 
luminosity of these clusters (iup), however, depends 
strongly on age. The light-to-mass ratio, or the flux 
of a cluster of constant mass, can be approximated as 
with C ^ 1.0/0.7/0.6 for the U/V/K-hsmd. Note 
that if A/up >> ^/max then the luminosity of clusters 
with Mmax (imax) wiU actually (statistically) increase 
with age (in log(age) bins) as Lmax oc t^~'' . The fact 
that the most luminous cluster in a galaxy is generally 
young therefore already argues for a truncation of the 
CIMF or a steepening at high masses. With Lup oc 
we know the change of Lup with age dLup/dt oc t~'^^^. 
Combining this with the constant diV(A/up)/d< we can 
thus write 



dN 
dL„n 



oc t' 



oc L 



C+i 



-1-1/C 

up 7 



(10) 

(11) 



-1/C 



where t is substituted for Lup ' 

We thus predict that the bright end of the LF is a 
power-law with index —1 — 1/0.7 = —2.5 in the V"-band, 
whereas the faint end is a power-law with index —2. 
This is in perfect a gree ment with what was fo und by 
Gieles et al.l (|2006bl ) and lWhitmore et al.l (|l999l ) for the 
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Fig. 3 Simulated cluster population with dN{M > 10^Mo)/df = lOMyr"^ and Mup = Mq. Top left: Age-mass 
diagram of all cluster brighter than My = —7. Top right: Age-A'Iv diagram. Bottom left: Mass distribution. Bottom right: 
Magnitude distribution. 



LF of clusters in M51 and the Antcnnacll], respectively. 
In the ?7-band ~ 1) the index of the LF is thus —2 
and a truncation of the CIMF is not observable in the 
[/-band LF. This filter dependent power-law index is 
also found for the bright end of the LF o f M51 clusters 
( Gieles et all l2006al lbl: iHaas et al.ll2008[ ). though with 



low significance. 

If Mmax is not sampled up to Mup at all ages, the 
index of the LF will be somewhere between —2.5 and 
—2. We illustrate this in Fig. [3] with a Monte Carlo 
simulation. We sample clusters from a CIMF with 
M,ni„ = IO'^A'/q and M^p ~ 10^ Af0, assuming a con- 
stant CFR between lO^yr and lO^yr and no disruption. 
We apply a detection limit at My = —7. In the top 
panels we show the masses and luminosities as a func- 
tion of age. For this CFR Mmax=Mup at 100 Myr, 
since M„jax = Mmin * dN{M > 10^ MQ)/dt * At = 
IG^ X 10 X 100 = 10^ Mq. In the top left panel we 
illustrate the statistically expected Afmax as a function 
of log(age) as a dashed line. The points follow the line 
with some clusters above the line, which is allowed since 
Mmax is not a hard limit. For ages older than 10® yr 

-A^max — ■ 

In the bottom panels we show the resulting mass 
and luminosity distributions together with power-law 
approximations for dN/dM and dN/dL. The mass 
function has a hard cut-off at high masses and is heav- 
ily affected by the detection limit at low masses. The 
LF, however, is still a power-law, but with a slightly 
smaller index (~ —2.3). The L F is in fact best ap - 
proached by a double power- law ICieles et all l2006a[ ) . 



^However, at this conference Brad Whitmore presented the LF 
of Antennae clusters derived from HST/ACS data and the break 
in the LF seems to have disappeared. 



but in the situation where Mmax is not sampled up to 
Mup at all ages, the LF becomes only slightly steeper 
than the underlying CIMF. Note that there are other 
effects that can make the LF even steep er, but this is 
only when M,„ax — Mup (jLarsenl 120021 ). We empha- 
sise again that without a truncation the integrated LF 
should be a power-law with index —2. 



5 Conclusion 

The fact that for M^ax < 10*^ Mq the relation be- 
tween Mmax and N{M > 10'^ M0) is similar in seven 
different galaxies argues for universality of the CIMF 
in that mass range. Whether the truncation mass of 
10^ Mq is universal is yet unclear. There are observa- 
tional hints that M.p in the Antennae is a factor of four 
higher than in M51 |Gieles et al.ll2006al|bl ). Moreover, 



known (e.g. 


Miller et al. 19971 Maraston et al. 2004: 


Bastian et al. 


2006; Reikuba et al. 20071: Wilson et al. 


20061). However, these obiects seem to be different 



from "regular" star clusters based on some fundamen- 
tal properties. Above ~ 10^ Mp there is a relation 



betw e en cluster radius and mass (jKissler-Patig et al 



200d iHasegan et all l2005h . while below - 10«Mq 
there in no such correlation, neither fo r young clusters 
( Larsenll2004 IScheepmaker et al ]|2007|l. nor for old clus- 



ters ( e.g. iDiorgoyski and MevlanI Il994 : McLaughlinl 
2OOOI: [Jordan et al.ll2005h . In a ddition, mass-to-light 



ratios increase above ~ 10^ Mq (jReikuba et al.l 120071 ) 
and the chem ical properties abo ve this mass are dis 



tinct as well ( Mieske et al. 20061) . Lastly, a universal 



upper mass helps to ex plain the universal g lobular clus- 
ter turn-over location (j Jordan et al.ll2007l ). 
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